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Preliminary results of a work in collaboration with M.R. Pennington are presented. 
Extending a scheme introduced by Tornqvist, we investigate a dynamical model in 
which the spectrum of scalar mesons can be derived, with the aim of locating the 
lightest glue-state. Adding hadronic interaction contributions to the bare propaga- 
tor, to 'dress' the bare quark-model qq states, we are able to write the amplitudes 
and the phase shifts in the approximation in which scalar resonances decay only 
into two pseudoscalar channels. The fit of these quantities to experimental data 
gives a satisfactory understanding of how hadronic interactions modify the under- 
lying 'bare' spectrum. In particular, we examine the case in which a glue-state is 
introduced into the model. 



1 Introduction 

Gluons carry colour charge which means they interact together. Consequently 
it is possible for them to cluster and form objects which are colourless overall. 
These would be just like conventional hadrons, but with constituents that are 
massless gauge bosons. These are known as glueballs. 

Interest in glueballs has increased since new candidates for gluonic states have 



emerged from experimental resultsllrEl especially in the 1.5 — 2.0 GeV energy 
region. Moreover, Lattice QCD calculations (in quenched approximation) have 
recently suggested the presence of light scalar glue-states in this same energy 
region 0. From the experimental point of view, it is now clear that there are 
too many confirmed scalar (0 ++ ) mesons to form one qq meson nonet. The 
problem is particularly pronounced for the 1=0 sector, since at least four /o's 
now appear in the particle data listing. As a consequence, we can infer that 
some of them have to be extra states. 

The quark model gives a reasonably good description of the vector and tensor 
meson spectra and properties, but its predictions for the scalar sector are very 
disappointing. To understand how and why scalars are so different from vec- 
tors and tensors, we consider a simple model in which all bare meson states 
belong to ideally mixed quark multiplets. We call rin the nonstrange light 
state and suppose that substituting a strange quark for a light one increases 
the mass of the state by Am s ~ 100 MeV, as illustrated in Fig. 1. 
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Figure 1: Spectrum of the bare qg states 
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The bare propagator for each of these bound states will be of the form I 

1 

Mf-s 

with a pole on the real axis, corresponding to a non decaying state, for example 

\4>)o = \ss) 

for the vector 1 = state, and 

l/o)o = |*5) 

for the scalar 1 = state. If we now assume that the experimentally observed 
hadrons are obtained from the bare states (nn, sn, ss, ...) by dressing them 
with hadronic interactions, the propagator becomes 

P(S) = M 2 (s)~s-iM(s)T(s) ' (2) 

and the pole moves into the complex s-plane. The corresponding state can be 
decomposed as 

\4>) = v 7 ! - e 2 |ss) + ei\KK) + e 2 \pn) + ... (3) 

where calculation would give e 2 = e 2 + e 2 , + ... -c 1. The hadronic loop con- 
tributions allow the bare states (ss in this example) to communicate with all 
hadronic channels permitted by quantum numbers, and this enables the 4> me- 
son to decay. A similar picture works for the tensors. 

For scalars the situation is different because the dominant decays are just into 
two pseudoscalars, the couplings are bigger and they couple strongly to more 
than one channel, creating overlapping and interfering resonance structures. 
Furthermore, being S-waves, the opening at thresholds produces a more dra- 
matic s-dependence in the propagator. As a consequence, the /o(980), for 
example, turns out to be predominantly a \KK) state, and not an |ss) one: 

|/o(980)> = Vl-e 2 \KK)+ei\ss) + ... . (4) 
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This does not mean that the /o(980) is a KK molecule, because the seeds 
of the model are conventional qq states and the binding forces are not due to 
inter-hadron interactions alone. The crucial point here is the fact that hadronic 
interactions allow the qq bare states to communicate with all possible hadronic 
channels and these channels to communicate with each other, giving rise to a 
mixing which means, for example, that the /o(980) spend most of its time in 
a KK state and not in an ss one. 



2 A closer look at the model 

Let us now examine the model in more detail starting, for simplicity, from the 
case in which just one resonance is produced. If we define a vacuum polariza- 
tion function n(s) which takes into account all the possible two pseudoscalar 
loopj-contributions to the propagator P(s), we can easily write its imaginary 
partta: 

Imn( S ) = = -^^^( S - S ^)^ 2 ( S )^( S - StM ) (5) 

i i 

where the index i runs over the pseudoscalar channels, the g^s are the SU(3) 
flavour couplings, the fej's are the cm. momenta and (s — sa,j) are the Adler 
zeros. F i (s) are the form factors, which take into account the fact that the 
interaction is not pointlike but has a spatial extension. These are parametrized 
by 

F < = exp (=W) • (6 > 

where the momentum fco is inversely proportional to the range of the interac- 
tion. The real part of the vacuum polarization function can be found from the 
dispersion relation 

Ren( S ) = IP r M ^1 (7) 
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No subtraction is needed, since the form factors decrease fast enough when 
\s\ — *■ oo. At this point we are able to write the propagator as 

1 

771*5 + n ( S ) 

and the contribution to the i —> j amplitude as 



P(») = ^2 , ^ I - (8) 



n ' ml + II(s) - s K ' 
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Notice that i?y(s), having its numerator and denominator related to each 
other by the same Gi couplings, is not the most general amplitude satisfying 
unitarity. This would look like: 

A i:j (s) = i?y(s)e 2lQ(s) +sina(s) e la{s) (10) 

where a(s) is an unknown function of s real along the right hand cut, and the 
second term in Eq. (8) accounts for the background contribution. To avoid 
an increasing number of parameters in the model, we will consider just the 
approximate amplitude Fig of Eq. (^), having checked that a constant value 
of a of about 15° in Eq. ([To]) would typically work just as well. 

Fig. 2 shows the behaviour of the real and imaginary part of the running 
complex mass function m 2 (s) = m 2 , + II(s), for the I = 1/2, A'q(1430), and 
the 1=1, ao (980). Here the parameters of the model (the mass of the bare 
rin state, the form factor cut-off ko, an overall coupling 7 and the position 
of the Adler zero) are determined by fitting the amplitudes and the phase 
shifts to the experimental data from LASSla. As we can see, the opening of 
each threshold gives an extra contribution (square root cusp) to the imaginary 
part of the vacuum polarization function. This is reflected in the shape of its 
real part, which presents a very strong s-dependence at every threshold. The 
intersection of Rem 2 (s) with the curve s represents the square of the Breit 
Wigner mass, and the value of Imll(s) at this point is related to the Breit 
Wigner width by 

r w = ±!!W. (ii) 

rriBW 

Notice how the negative contribution of Reil(s) shifts down the actual mass 
of the resonance with respect to the value of the corresponding bare mass. 
This effect is particularly pronounced for the ao(980) because the mass of the 
dressed bound state happens to coincide with the first of two thresholds, the 
A A and the 7777' ones, which are very close to each other and have similar cou- 
plings: despite the bare mass mo being fixed at 1420 MeV, the Breit Wigner 
mass of the oq is found to be as low as 987 MeV. 
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Figure 2: Real and imaginary parts of the complex running mass functions for the (a) 
^(1430), and (b) a (980), illustrating the effect of thresholds. 



To analyze the more controversial I — scalar sector, we need to examine the 
case in which more than one resonance is produced. Now, the polarization 
function has to incorporate the features (couplings, masses, form factors, ...) 
of each of the intermediate particles created in the process: 

Imn Q(3 ( S ) = -Y,G ati {s)GiA*) - ( 12 ) 

i 

ReU a ,(s) = l -P r ds'^^l, (13) 

77 Jsth,i S S 

where the indices a and /3 run over the N different resonances. As a conse- 
quence, the propagator 

Pafiis) = -1-^ rrr- tt 7"T ( 14 ) 

11 {ml- s)(5 Q/3 +n a/3 (s) 

becomes an N x N complex matrix, and the couplings are iV-dimensional col- 
umn vectors. The amplitudes are determined using a diagonalization procedure 
and can be written as 



~ m a,diag( S ) 



where the diagonalized masses m a ,diag(s) are admixtures of all the bare masses 
and the vacuum polarization matrix elements. They are, like the new couplings 
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G' ai (s), complex and s-dependent. 

What is the interpretation of this more complicated picture ? The physical 
observed hadrons are the states we obtain after diagonalization, an admixture 
of the seed states of the model. The mixing among them is embodied in the 
diagonalization procedure, which allow all the channels to communicate with 
each other and to create new physical states with masses, couplings and widths 
different from the ones of the primitive states. 



We now want to consider explicitly the case in which not only the two conven- 
tional fo and /q are present in the I — sector, but also a third state arises, 
thanks to the presence of a glue gg seed. Since, in principle, we do not know 
the mass of the bare glue-state, m gg is assumed to be an extra parameter of 
the model. To avoid a further increase in the number of parameters we will 
also assume that the gg state cannot mix with the qq states at the bare level: 
all the mixing will occur via hadron interactions. The bare couplings for the 
glueball to the two pseudoscalar channels are the ones of an SU(3) flavour 
singlet. Fig. 3 presents the results we obtain if we choose a gg bare mass of 
1.8 .QeV and readjust the other parameters to fit the data from irir scatter^ 
ingEEliil, in such a way that the agreement with the LASS data on Kq (1439)0 
remains satisfactory. The two lower states have the_same features that would 
be obtained if only two resonances were consideredu: a very broad state with 
a mass of about 1 GeV and a narrow one with a mass of about 1.2 GeV. 
The higher glue-state almost decouples from the lowest channels, but shows 
relatively strong couplings not only to the heaviest rj'rj', but also to the r\rj 
channel, to which it had a zero underlying bare coupling. Again, the Breit 
Wigner mass is considerably lower then the input one: for these values of the 
parameters our glueball mass is about 1.57 GeV and its width is about 200 
MeV. Notice how the masses of all the scalar 1 = mesons are related to the 
positions of the two pseudoscalar thresholds. The amplitude presents two dips: 
the former corresponds to the presence of a very narrow resonance sitting on 
top of a very broad one, the second is related to the third heavier resonance, 
in agreement with the experimental data from Crystal Ball El, GAMSE3, and 
the analysis from Bugg et al. EJ and Anisovich et al. Il3 
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Figure 3: (a) Real parts of the diagonalized masses of the three / = scalar resonances 
created in tttt scattering, (b) Modulus of the amplitude corresponding to the same 
process 



3 Conclusions 

In this talk I have presented some preliminary results of a work performed in 
collaboration with M.R.Pennington, which is still in progress. In particular, a 
proper fit of the more recent experimental data up to about 1.8 GeVE3E3 has 
still to be completed. 

The model we are using requires approximations, to allow the calculation to 
be performed. The amplitudes are assumed to be pole dominated and so the 
militarization chosen is not the most general one. Furthermore, only decays 
into two pseudoscalars are considered, whereas other thresholds like multip- 
ion, vector-vector or axial-pseudoscalar ones are neglected. 
Despite the assumptions, this scheme has many great advantages: first of all 
it gives a dynamical description of the particularly complex mechanism which 
leads to the creation of scalar resonances, naturally taking into account the 
mixing amongst different states, and explaining why the scalars differ from 
vector and tensor mesons. Moreover, it is simple and with very few parame- 
ters. 

Consequently, we believe this model is not a machine from which one can 
blindly extract numbers from a fitting program, but an efficient schematic way 
to approach the dynamics of the non-perturbative hadronic world: a world in 
which so much remains to be understood. 
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